Quantum- like approaches to the beam halo 

problem * 



Sameen Ahmed KHAN 
Dipartimento di Fisica Galileo Galilei Universita di Padova 
Istituto Nazionale di Fisica Nucleare (INFN) Sezione di Padova 
Via Marzolo 8 Padova 35131 ITALY 



E-mail: khan@pd.infn.it, \http: / /www. pd.infn.it/~khan/ 



Modesto PUSTERLA 
Dipartimento di Fisica Galileo Galilei Universita di Padova 
Istituto Nazionale di Fisica Nucleare (INFN) Sezione di Padova 
Via Marzolo 8 Padova 35131 ITALY 



E-mail: pusterla@pd.infn.it, http://www.pd.infn.it/~pusterla/ 



Abstract 

An interpretation of the the "halo problem" in accelerators based on 
quantum-like diffraction is given. Comparison between this approach and 
the others based on classical mechanics equations is discussed 
Keywords: Beam Physics, Quantum-like, Beam halo, Beam Losses, Stochasticity. 



I. INTRODUCTION 

Recently the description of the dynamical evolution of high density beams by using the 
collective models, has become more and more popular. A way of developing this point of 
view is the quantum-like approach [0] where one considers a time-dependent Schrodinger 
equation, in both the usual linear and the less usual nonlinear forms, as a fluid equation 
for the whole beam. In this case the squared modulus of the wave function (named beam 
wave function) gives the distribution function of the particles in space at a certain time @ . 
The Schrodinger equation may be taken in one or more spacial dimensions according to the 
particular physical problem; furthermore the fluid becomes a Madelung fluid if one choses 
the equation in its usual linear version. 



* To appear in: Proceedings of the 6th International Conference on Squeezed States and Uncer- 
tainty Relations, (ICSSUR'99), 24-29 May 1999, Napoli, ITALY, (NASA Conference Publication 
Series). 
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Although the validity of the model relies only on experiments and in particular on new 
predictions which must be verified experimentally, we like to invoke here a theoretical argu- 
ment that could justify the Schrodinger quantum-like approach we are going to apply. Let 
us think of particles in motion within a bunch in such a way that the single particle moves 
under an average force field due to the presence of others and collides with the neighbouring 
ones in a complicated manner. It is obviously impossible to follow and predict all the forces 
deterministically. We then face a situation where the classical motion determined by the 
force-field is perturbed continuously by a random term, and we have a connection with a 
stochastic process. If one simply assumes that the process is Markovian and Brownian, one 
obtains following Nelson ||, a modification of the classical equations of motion that can 
be synthesized by a linear Schrodinger equation which depends on a physical parameter 
having the dimension of action 0. Usual wave quantum mechanics follows if this parame- 
ter is chosen as the Planck's constant H, whereas the quantum-like theory of beams in the 
linearized version is obtained if one choses the normalized emittance e 0. In both cases, 
quantum particle and beam respectively, the evolution of the system is expressed in terms 
of a continuous field ip defining the so called Madelung fluid. We may notice that the nor- 
malized emittance e having the dimension of action is the natural choice for the parameter 
in the quantum-like theory that corresponds to Planck's constant h in the quantum theory, 
because it reproduces the corresponding area in the phase-space of the particle. 

We here point out that, after linearizing the Schrodinger-like equation, for beams in an 
accelerator one can use the whole apparatus of quantum mechanics, with a new interpre- 
tation of the basic parameters (for instance the Planck's constant % — > e where e is the 
normalized beam emittance) and introduce the propagator K (xf,tf\xi,ti) of the Feynman 
theory for both longitudinal and transversal motion. A procedure of this sort seems partic- 
ularly effective for a global description of several phenomena such as intrabeam scattering, 
space-charge, particle focusing, that cannot be treated easily in detail by "classical mechan- 
ics" and are considered to be the main cause for the creation of the Halo around the beam 
line with consequent losses of particles. 

Let us indeed consider the Schrodinger like equation for the beam wave function 

iedt^ = ~^ + U(x,t)^ (1) 
2m 

in the linearized case U (x,t) does not depend on the density \ipf- e here is the normalized 
transversal beam emittance defined as follows: 

e = m cype, (2) 

e being the emittance usually considered, where as we may introduce the analog of the de 
Broglie wavelength as A = e/p. We now focus our attention on the one dimensional transver- 
sal motion along the x-axis of the beam particles belonging to a single bunch and assume a 
Gaussian transversal profile for a particles injected in to a circular machine. We describe all 
the interactions mentioned above, that cannot be treated in detail, as diffraction effects by 
a phenomenological boundary defined by a slit, in each segment of the particle trajectory. 
This condition should be applied to both beam wave function and its corresponding beam 
propagator K. The result of such a procedure is a multiple integral that determines the 



2 



actual propagator between the initial and final states in terms of the space-time intervals 
due to the intermediate segments. 

K (x + x Q ,T + r\x', 0) = / K (x + x , t\x + y n ,T + (n - l)r') 

J-b 

xK(x + y n ,T+(n~ l)r'\x + y n _ u T + (n - 2)r') 

x • • • K (x + yx, T\x', 0) dy x dy 2 ••■dy n (3) 

where r = nr' is the total time of revolutions T is the time necessary to insert the bunch 
(practically the time between two successive bunches) and (—6, +b) the space interval defin- 
ing the boundary conditions. Obviously b and T are phenomenological parameters which 
vary from a machine to another and must also be correlated with the geometry of the vacuum 
tube where the particles circulate. 

At this point we may consider two possible approximations for K (n\n — 1) = 
K 0&o + Vn, T + (n - 1)t'\x q + y n „x + {n- 2)r'): 

1. We substitute it with the free particle K assuming that in the r' interval (V r) 
the motion is practically a free particle motion between the boundaries (—6, +b). 

2. We substitute it with the harmonic oscillator (n\n — 1) considering the harmonic 
motion of the betatronic oscillations with frequency uj/2h 



II. FREE PARTICLE CASE 



We may notice that the convolution property (|3|) of the Feynman propagator allows us 
to substitute the multiple integral (that becomes a functional integral for n — > oo and 



t' — > 0) with the single integral 

h6 



K (x + x , T + t\x', 0) = ( dyK(x + x ,T + t\x + y,T) K {x + y,T\x' ,0) dy (4) 

J — b 

In this note we mainly discuss the case 1. and obtain from equation (|J) after introducing 



the Gaussian slit exp 
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instead of the segment (—6, +b) we obtain from 
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where vo = and xois the initial central point of the beam at injection and can be 

chosen as the origin (a;o = 0) of the transverse motion of the reference trajectory in the test 
particle reference frame. Where as h must be interpreted as the normalized beam 
emittance in the quantum-like approach. 

With an initial Gaussian profile (at t = 0), the beam wave function (normalized to 1) is 



3 



/(*) 



1 

a l 4 



7T 



cxp 



a 



12 



r.m.s of the transverse beam and the final beam wave function is 

1 

f + OO 



[x 



dx'(-Y ei-^K (x, T + t; x', 0) = B exp 

oo \7T/ 



(6) 
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The final local distribution of the beam that undergoes the diffraction is therefore 
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p(x) = \4>(x)\ = BB*exp 
where a = — (C + C*) and the total probability per particle is given by 

r+oo pff 

P = / dxp(x) = BB*J- 
J-oo V a 
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One may notice that the probability P has the same order of magnitude of the one computed 
in || if ^= is of the order of b. 



III. OSCILLATOR CASE 



Similarly we may consider the harmonic oscillator case (betatronic oscillations) to com- 
pute the diffraction probability of the single particle from the beam wave function and 
evaluate the probability of beam losses per particle. The propagator K w (x,T + r\y,T) in 
the later case is: 
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where 
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(14) 



p u (x) = \<PM\ 2 = N* N exp - (M* + M)x 
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From the approximate formulae [Tl] and |T^ we notice that the parameter r does not play 
a significant role in the calculation of the probabilities. We gave it a value r = 1 sec, 
considering about 10 6 revolutions in LHC and HIDIF storage rings. 



IV. PRELIMINARY ESTIMATES 



Parameters 

Normalized Transverse Emittance 
Total Energy, E 
T 
b 
P 



TABLE-I: Free Particle Case 
LHC 

3.75 mm mrad 
450 GeV 
25 nano sec. 
1.2 mm 
3.39 x 1(T 5 



HIDIF 

13.5 mm mrad 
5 Gev 

100 nano sec. 
1.0 mm 
2.37 x 10~ 3 



TABLE-II: 

Parameters 

Normalized Transverse Emittance 

Total Energy, E 

T 

b 

u 



Oscillator Case 
LHC 

3.75 mm mrad 
450 GeV 
25 nano sec. 
1.2 mm 
4.47 x 10 6 Hz 
3.44 x 10~ 5 



HIDIF 

13.5 mm mrad 
5 Gev 

100 nano sec. 
1.0 mm 
1.12 x 10 7 Hz 
2.96 x 10~ 3 
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V. CONCLUSION 

The parameters entering into the probability formulae are very few and must be looked 
at as purely phenomenological. To be more specific, b, r, and T (b in particular) hide 
several fundamental processes that may be present in the beam bunches and that play a 
deterministic role in the creation of the halo such as intrabeam scattering beamstrahlung, 
space-charge effects and imperfections in the magnets of the lattice that cause nonlinear 
perturbative effects. 

The fact that such a small amount parameters take into account many physical processes 
is a nice feature of the quantum-like diffraction approach. However a way of connecting this 
method with the physical processes mentioned above as well as with the nonlinear dynamical 
classical theory is mandatory at this point. 

Another interesting feature of the parameters used is that their numerical values are very 
reasonable because they are within the ranges. One expects: T may be considered as the 
average time interval between the two successive injection, r the time interval between two 
successive diffractions (r = tit 1 is the total time of revolutions) and 2b the phenomenological 
diffraction slit width. We recall that in the usual optics diffraction through a slit is also a 
macroscopic means of dealing with many complicated physical effects such as scattering of 
light, electrons etc., at the atomic level. 

The two relevant concluding remarks are the following: 

1. The probability calculated from the free and the harmonic oscillator propagators (both 
in the transversal motion of the particles) appear very close for the two different circular 
systems such as LHC and HIDIF rings. 

2. The HIDIF scenario, as expected has a total loss of beam power which is at least 10 3 
times higher than LHC. 

These preliminary numerical results are encouraging because they predict halo losses 
which seem under control. Indeed the HIDIF scenario gives a total loss of beam power per 
meter which is about a thousand higher than the LHC; however in both cases the estimated 
losses appear much smaller than the permissible 1 Watt/m. 
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